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FINITE ELEMENT MODELS OF WIRE ROPE 


FOR VIBRATION ANALYSIS 

J. E. Cochran, Jr.,* N. G. Fitz-Coy^ and M. A. Cutchins* 

Auburn University, Alabama 

Abstract 

The usefulness of wire rope in shock and vibration isolation is briefly 
reviewed and its modeling for the purpose of vibration analysis is addressed. 

A model of a nominally straight segment of wire rope is described in which the 
rope structure is represented by a maiden, or central, strand of wire with one 
(or more) strand(s) wrapped around it in a helix (helices). The individual 
strands are modeled using finite elements and MSC NASTRAN. Small linear 
segments of each wire are modeled mathematically by dividing them lengthwise 
into triangular prisms (thick "pieces of pie") representing each prism by a 
solid NASTRAN element. To model pretensioning and allow for extraction of 
internal force information from the NASTRAN model, the "wound" strands are 
connected to the maiden strand and each other using "spring" (scalar elastic) 
elements. Mode shapes for a length of wire rope with one end fixed to a 
moving base and the other attached to a "point" mass, are presented. The use 
of the NASTRAN derived mode shapes to approximate internal normal forces in 
equations of motion for vibration analyses is considered. 


^Professor, Aerospace Engineering 
tGraduate Research Assistant 
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Introduction 


Wire rope 1 is, from the basic point of view, simply several strands of 
wire twisted, or wound, together (see Fig. 1). Some types are commonly called 
"cable" and are used to carry electricity, support bridges and "cable cars," 
raise and lower heavy loads and in many other practical ways. A less obvious, 
but equally important, use of wire rope is in shock and vibration isolation 
devices. 2 » 3 The structure of wire rope provides many interfaces at which a 
portion of the relative motion of strands of wire is converted by friction 
into heat, thereby dissipating vibrational energy. 4 Furthermore, the 
stiffness of wire rope structures can be tailored to provide support and 
restoring forces. Stiffness and damping are adjusted by varying wire 
diameter, the number of strands, pretensioning and the arrangement of lengths 
of the wire rope. Commonly, helical coils of ropes 2 are fixed in clamps 
(see Fig. 2) to form individual shock and/or vibration isolators. The 
isolators are used to support and isolate communications equipment in vehicles 
which are subjected to large magnitude, short-term accelerations; i.e., 
"shocks." In addition to absorbing shock, the internal, or system, 
damping, *♦» 5 of the wire rope devices provides vibrational isolation over wide 
ranges of frequencies and amplitudes. 

The damping characteristics of wire rope and vibration isolators made 
from it are not well understood from the theoretical standpoint. Apparently, 
the design of individual isolators is accomplished by experimentation by 
engineers with considerable experience in applications of these devices. 2 
Realistic mathematical models of wire rope isolators would be useful in the 
design process and perhaps would allow the achievement of the confidence 
levels in isolator characteristics needed for more applications in which 
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Fig. 1. Wire rope. 


1254 




damping rates and dynamic response must be very accurately known to prevent 
resonance and control interaction problems. 

The purposes of this paper are (1) to present the results of work on the 
problem of obtaining accurate mathematical models of wire rope isolators using 
MSC NASTRAN and (2) to present some preliminary results of work toward 
obtaining equations of motion of a wire rope "pendulum" for use in correlating 
experimental and theoretical results. 

Wire Rope Model 

Basic Model 

A sketch of a wire rope” segment of length L, composed of a maiden wire 
strand and a single "wound” strand, is shown in Fig. 3. It is assumed that 
when the rope is in its undeformed state the angle of the helix formed by the 
centerline of the wound strand, a, is constant over the length L. For the 
finite element analysis, the maiden strand is divided into 372 elements 
such that each cross section of the strand looks like a hexagon cut into tri- 
angular pieces which have thicknesses equal to the radius of the strand. The 
other strands are modeled in a similar manner. In addition, scalar elastic 
elements ("springs”) are used to connect the wound and maiden strands. These 
springs are incorporated so that the distribution of the normal force between 
the two strands can be determined for each mode shape computed (see Fig. 5). 
NASTRAN produced figures showing two-strand and seven-strand wire rope are 
given as Figures 6 and 7, respectively. 

Boundary Conditions 

The boundary conditions are chosen so that the "top" of the segment is 
fixed. Two choices for the boundary conditions on the other end are 
considered here. The first is to use a "free" end condition, (Case 1). The 
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Fig. 7. NASTRAN model of seven-strand wire rope. 
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second choice (Case 2) is to attach to the "bottom" end a mass which is 
relatively large compared to the mass of the rope segment. 

The fourth mode shape for the seven-strand rope alone (Case 1) is shown 
in Fig. 8. The corresponding mode shape for the Case 2 "pendulum" model is 
given in Fig. 9. 


Equations of Motion 

Our stated objective is to obtain "analytical" models of wire rope 
vibration isolation. Here, "analytical" implies that mode shapes from NASTRAN 
models without damping will be used rather than analytical mode shapes such as 
the classical ones for beams. The mode shapes are to be used to approximate 
the displacements of the wire rope and attach "rigid" bodies when linear 
damping is present and when nonlinear nonconservative and conservative forces 
act on the system. The wire rope "pendulum" is a simple system which we 
hope to use to test this method for modeling these effects. 

Except for the determination of nonconservative internal forces, the most 
attractive way of deriving equations of motion for the pendulum is to use the 
Lagrangian formulation. 

The Lagrangian 

To obtain the Lagrangian ^ , we must write the kinetic and potential 
energies of the pendulum. The support point. A, is assumed to move in a 
prescribed manner so that the vector R A from the fixed point 0 is known. The 
kinetic energy, T, of the pendulum, assuming the mass M is a "point mass," is 

"I,", / < i* + £d» ) -<iA + id. )d "‘ 

i-1 

+ 1/2 M (V A + r M ).(V A + r M ) (1) 
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Fig. 8 Mode shape for wire rope model. Case 1. 
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Fig. 9 Mode shape for wire rope model 


Case 2. 


where V^* , the velocity of A, is the velocity of a generic mass 

element dm, N is the number of strands (see Fig. 3), and m* is the mass of the 

m 

fth strand. We approximate r, by £ j>. q., where the j>. are mode shapes, 

“ m j=l J 3 J 

the qj are generalized coordinates, m is the number of modes used and 
m 

r w * Y d>. (0 0 L)q,. The A. are functions of the coordinates x , y and z 
-M jtl j O* J o o 

of dm in the undeformed "structure" and the qj are, of course, functions of 
time to be determined. 

In matrix form, Eq. (1) is 

T - ? “l V + M si k i + 7 i T &3 M 


+ E / * £ d® + I 7 / £ T | T | £ dm 


:a l i s a i 2 

i“l m^ i*i m^ 


where Mf is the total mass of the system, 3 = (q^q 2 * • •£ m ) T 

and $ = an d * 8 = evaluated at x q = y Q = 0 , z q ■ L. 


( 2 ) 


As expected, we may write 

„ 1 »T „ • , T • , , 

T = j j M 3 + a 3 + b , 


(3) 


where 


N T T 

M = l J | | dm + M | L | l 


i=l m. 


(4a) 


and 


a T = MV*# T + V? I I i dm 
- -A=L -A . L , J =■ 
i=l m^ 


» -Hv 


(4b) 


(4c) 
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the potential energy, V, is due to the elastic stiffness of the system and to 
gravity. For the former, write 

v s - 1/2 a T g s S (5) 

where K is the stiffness matrix. 

=s 

The potential energy due to the location elements of mass in the rope and 
the mass M in a uniform gravity field is obtained by assuming that M is a 
point mass at the end of the length of rope and that the rope has mass per 
unit length c(z 0 ). Furthermore, the motion is restricted to the xz-plane and 
dx/dz 0 = (100) (d£/dz 0 <j) is assumed small enough for terms of order q3 and 
higher to be neglected. Then, the potential energy contribution due to 
gravity can first be written as 
L 

Vg ” / Mg(l - cos 0)dz Q 
L z 0 

+ / o[l - / cos[0(5)]d?]dz (6) 

o 0 

By noting that cos0 s 1 - y (dx/dz Q ) 2 , we may write 



Obviously, V a is of the form, 



The nonconservative forces acting on dm are assumed to be due to 
viscosity and coulomb friction between the strands. Viscous forces can be 
included in the equations by using some form of modal damping. For the coulomb 
friction, we need to consider the rope's structure. 

Let the strands be renumbered so that the maiden strand is i=0, and at 
any value of z Q , the other strands are numbered 1 through N-l in a 
counterclockwise manner about the z 0 axis. Also, let f^ denote the force on 
denote the force on the element dm^ of the i fc h strand and let p denote the 
coefficient of friction. Furthermore, let fjj denote the magnitude of the 
normal force between the i th and the j th strands. Then, for dmo, we have the 
coulomb friction force, 

■ A u - W- ( Mi + <U) 

O 1*1 0 1 


where t^ and bj are the tangent and bi-normal unit vectors, respectively, of 
dm^. For dmj, j > 0, we may write. 


-dm. 


+ “ f oj as “[ ( Ed. o - + ijty 


- uf Jk •nI<E d . j - Ed^-'EjEj + SjV 


•j — t 

with j, k and i in the following triplets: 


j 


— " 
1 


Y 


’3' 


Y 


N-l 

k 

■- 

2 

9 

3 

> 

4 

9 

5 

9 • • * 9 

1 



6 


1 


2 


3 


N-2 


( 12 ) 


(13) 


In Eqs. (12) we have also used** 
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(14a) 


tj - cos a [- sin 0^ i + cos 0jjJ + sin a k 

and 


- b j 


A A A 

sin a(sin 0^1- cos O^j] + cos a_k 


(14b) 


where 0j(z o ) is the angle of the centerline of the j 1 -* 1 strand at z Q . 

Because r dm - R A + |(x o ,y o ,z Q ) a» where (x o ,y o ,z Q ) is a point in the 
undeformed i^* 1 strand, dr^ * £ (x o »y o »z o )6a- 8 enera li ze d force matrix is 

therefore 7 


N-l _ 

2-1 / a T 

i-0 V 

where the integral indicates a summation of the quantities within it. 
The equations of motion then follow from « T-V s -Vg, viz.. 


(15) 


dt 


^ - 12 ?, 


3a 


3£ 


2 t - (2 a> T » 


(16) 


where - D g is the generalized viscous damping force. For normal modes of the 
undamped structure, we have (with D a diagonal matrix). 


Mj qj + kjqj + K^a - Qj - - k j (l 7) 

where Mj is the generalized mass and kj the generalized stiffness, 
respectively, for the mode is the row of K^, is the jth element 
°f 2 , d j tl »e diagonal element of D in the row and a^ is the j th element 


The motion in qj is coupled in two ways. First, through the gravity 
terms and, second, through the Qj since up to m of the appear in to 
determine its sign. 

The next steps to be taken are (1) to obtain experimental data for the 
position of M and (2) to match it with 

5m " h * Sm <l8) 

which can be determined when is specified and the q^ are known. 

Summary 

The use of wire rope in shock and vibration isolation devices has been 
reviewed briefly. Finite element models of a nominally straight length of 
wire rope have been described. These models were developed as part of an 
effort to construct mathematical models of wire rope for use in the analysis 
of vibration isolation devices which are constructed from wire rope. 

Equations of motion for one of the finite element models, a seven-strand 
rope suspended from one end and with a mass attached to the other end to form 
a "pendulum," were derived. The use of these equations to simulate the motion 
of the pendulum was discussed. 

Further work needs to be done with the current models and new models of 
helical isolators should be developed. 
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